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Abstract. The conversion of the read-out from the anti-symmetric port of the 
LIGO interferometers into gravitational strain has thus far been performed in 
the frequency domain. Here we describe a conversion in the time domain which 
is based on the method developed by GEO. We illustrate the method using the 
Hanford 4km interferometer during the second LIGO science run (S2). 
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1. Introduction 

The LIGO interferometers Ej are part of a world-wide network of gravitational 
wave detectors constructed to detect waves from astrophysical sources such as spinning 
neutron stars, coalescing binary black holes and neutron stars, and supernovae. 

Each of the two arms of the LIGO interferometers forms a resonant Fabry- 
Perot cavity. This increases the effective arm-length and thus the sensitivity of the 
instruments. Gravitational waves (along with seismic and other noise) change the 
relative length of the optical cavities and produce an external strain, 



that is incident on the interferometer. Here, L x {t) and L y (t) are the effective lengths 
of the x and y-arms, respectively and Lq the effective length of the cavities in the 
absence of an external strain. 

The output of the interferometers, however, is not h(t). Rather, the output is 
derived from light that escapes the anti-symmetric ("dark") port and provides the 
gravitational wave read-out. It is often called the error signal or "gravitational wave 
channel". Here we will refer to it as q(t). 

The re-construction of the gravitational wave strain hit) from the error signal 
q(t) is an essential part of the data analysis. Since the LIGO instruments are linear 
the gravitational wave strain incident on the interferometer is a a linear functional of 
q(t), namely, 



where R(t) is a suitable convolution kernel. Convolution in the time domain is simple 
multiplication in the frequency domain so that 



where x(f) denotes the Fourier transform of x(t). The ratio of the gravitational 
wave strain to the error signal is called the response function R(f). Finding the 
response function, or alternatively the kernel R(t), is the procedure referred to as the 
calibration. 

So far LIGO has used a calibration implemented entirely in the frequency domain 
@|, i.e. using Eq. On the other hand, the British-German interferometric detector 
GEO600 |S] has always produced a real-time calibrated h(t) using time-domain filters 
UJ- Some progress toward the production of h(t) for the LIGO interferometers has 
already taken place [7j. In this work we adapt the GEO time-domain calibration 
method to the case of the LIGO interferometers. 

In Section [21 we provide a simple description of the length sensing and control 
system of the LIGO interferometers, review the frequency domain calibration 
procedure and formally show how the strain can be reconstructed from the error signal 
of the interferometer in the time domain. In Section [3] we describe a method to track 
changes in the response using sinusoidal excitations injected into the instrument. In 
Section 21 we explain how the digital filters used to reproduce and invert the responses 
of various parts of the length control system are created and show their performance 
in the case of the Hanford 4km interferometer during the second science run (S2). We 
also describe our signal processing pipeline. We conclude in Section |SJ 
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Figure 1. Simple model for the feedback loop of the LIGO length sensing and 
control system. The three filters that make up the feedback loop are the sensing 
function C = a(t)Co, which converts the residual strain x r into the error signal 
q, the feedback filter D = j3(t)Do, which converts the error signal into the digital 
control signal d, and the actuation A which converts the digital control signal into 
the control strain x c that is subtracted from the external strain h. The digital 
signals q(t) and d(t) are recorded by the data acquisition system at 16384 Hz. 



2. Simple description of the LIGO length sensing and control system 

At low frequencies the external strain h is dominated by seismic noise and a control 
strain, x c , is subtracted from it (see Figure^ by physically moving the mirrors in the 
Fabry-Perot cavity to compensate for the seismic noise. This ensures that the residual 
strain that enters the optical cavity, x r , remains small at low frequencies and keeps 
the optical cavities in resonance. 

The optical cavity of the feedback loop is represented by the so-called sensing 
function, C = a(t)C'o. It converts the residual strain x r into the digital error signal q 1 
also called the gravitational wave channel, which is sampled at 16384 Hz and recorded 
by the data acquisition system. The error signal is measured in arbitrary units called 
counts. Here, Co is a reference sensing function that is measured at some time and 
o,{t) is an overall (real) gain that depends on the light power stored in the Fabry- 
Perot cavity which changes with the alignment of the mirrors. It is important to keep 
track of this gain and the method we use is explained in Section |3J The frequency 
dependence of the sensing function is determined primarily by the Fabry-Perot cavity 
in each arm and corresponds to a real pole at around 90 Hz. 

A digital feedback filter, D — /3(t)Do, is applied to the error signal q which 
produces a digital control signal d. The digital control signal, like the residual signal, 
is measured in units of counts. It is recorded by the data acquisition system at 
16384 Hz. Here, Dq is the feedback filter at some reference time and (3(t) is a real 
overall gain. The gain f3(t) is also recorded by the data acquisition system and may 
vary in time. For the first and second science runs (SI and S2), however, it was kept 
fixed. The filter Dq consists of a double pole at DC and various low pass filters as 
well as additional filters that keep the optics aligned and the Fabry-Perot cavity in 
resonance. 

The control signal d is converted to strain and used to adjust the length of the 
cavities via the actuation function A. This is achieved by converting the digital signal 
d into currents that run through coils placed near magnets glued on the back of the 
mirrors in the Fabry-Perot cavity. The frequency response of the actuation function 
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is largely determined by the pendulum suspension of the mirrors and thus consists of 
a pair of complex poles at the pendulum frequency f p w 0.75Hz. 

The re-construction of the external strain from the output of the interferometer 
feedback loop is referred to as the calibration. Calibration of LIGO data has thus far 
been performed in the frequency domain by constructing a response function R(f) 
that acts on the error signal q(f) of the feedback loop, namely, 

Hf)=R(f)q(f). (4) 

The response function can be derived from the feedback loop equations 

xr(f) = h(f)-x c (f), (5) 

q(f) = a(t)C (f)x r {f), (6) 

x c (f) = A(f)d(f), (7) 

d(f) = 0{t)D o (f)q(f). (8) 

From these equations it is easy to show that R(f) in Eq. 0} is given by 

l + a(t)(3(t)G (f) 
a{t)C Q {f) 

where Go = ACqDq is the reference open loop gain. The optical gain a(t) typically 
varies on a time-scale of seconds. The digital gain (3(t) may also vary but was kept 
fixed during SI and S2. Hence Eqs. H5I8|) apply for frequencies above a few tens of Hz. 

We would like to perform this procedure entirely in the time domain. To do this 
we reconstruct the strain time-series h(t) from the residual and control strains x r (t) 
and x c (t) as follows. The residual strain x r is given by 

x r (t) = h(t)-x c (t), (10) 

so that 

h(t) = x c (t) + x r (t). (11) 

The low frequency part of the external strain is dominated by the control strain x c (t) 
and the high frequency part is dominated by the residual strain x r (t). 

The optical gain of the sensing function a(t) typically varies on much longer 
time-scales than the decay time of the impulse response of the inverse of the sensing 
function. Thus, the residual strain can be formally re-constructed according to 

1 

a~ty 

where IU-i is a linear time-invariant filter operator for the inverse of the sensing 
function Cq. 

Similarly the control strain x c (t) can be constructed from d(t) by filtering it 
through the actuation function 

x c (t) = T A [d(t)}. (13) 

Therefore we can write the time series for the external strain in terms of the error and 
digital control signals, q(t) and d(t) respectively, as 

h(t)n^T c -r[q(t)]+T A [d(t)]- (14) 



*r(*)«— T c -i[q(t)}, (12) 



a(t) 



Although d(t) is an output of the interferometer servo loop recorded by the data 
acquisition system it can be computed from the error signal q(t) again assuming that 
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(3(t) varies in time much more slowly than the decay time of the impulse response 
function of the feedback filterj. In particular, 

d(t)^f3(t)T Da [q(t)}. (15) 
So in terms of just the error signal q(t) the external strain can also be written as 

Kt) w ^T c -i [<?(*)] + T A |)9(t) T Do [<?(*)]] . (16) 

Note that, because (3(t) is constant we could have factored it out of the second term 
on the right hand side of Eq. (jTB)l . 

Thus, in order to produce the time-series for the external strain we need 1) to 
track changes in the optical gain of the cavity, and 2) construct time-domain digital 
filters for the inverse of the sensing function, the actuation function and the feedback 
filter. These problems are dealt with in the next two Sections. 



3. Tracking changes in the calibration 

In order to track changes in the optical gain of the cavity, calibration lines are injected 
into the instrument by adding sinusoidal excitations to the control signal. This is a 
standard technique used in gravitational wave detectors. In the LIGO interferometers 
three calibration lines are injected; one at around 1kHz, another one around 150 Hz 
near the unity gain frequency of the open loop transfer function, and a third one at a 
few tens of Hz. 

The excitations may be added directly to the digital control signal d. Alternatively 
they may injected into one of the two arms, as shown in Figure[21 During S2 they were 
injected into the :r-arm of the interferometers and for the purposes of the following 
calculation we will assume this to be the case. Generalisations to injections into the 
y-arm or to the case where they are directly added to the digital control signal d are 
trivial. 

Figure |2] depicts a more detailed model of the actuation function A that shows 
the injection point. The actuation function A is given by 

A = K y g y A y - K x g x A Xl (17) 

where g x and g y are analog gains measured at DC that depend on the configuration of 
each of the arms and K x and K y are digital gains that are mainly used to compensate 
for differences in the analog gains. The functions A x and A y are the same for both 
arms except for digital notch filters that take out frequency components of the signal 
at the violin mode frequencies of the mirror suspensions. Since the calibration lines 
are not injected at those frequencies, for our purposes the functions A x and A y are 
identical along both arms and we can write 

A x == A y — A*. (1^) 

The excitation I is added to the x-arm signal after the digital gain K x is applied 
to the digital control signal. Therefore rather than Eq. Q, the feedback equation for 
the control signal at the frequency of a calibration line f c reads 

x c (f c ) = A(f c )d(f c ) - g x A4f c )l(f c ). (19) 

\ It should be noted that for the digital filter Do as it is currently defined this condition is never 
satisfied: The double pole at DC results in an ill-behaved digital filter that does not have a well- 
defined impulse response time. Later we show how to construct a modified digital filter D' which 
does satisfy this condition. 
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Figure 2. Details of the actuation and calibration line injection point for the 
LIGO interferometers length sensing and control system. This is the content of 
the actuation function A in Figure HI 



Using Eqs. (|17(1 and l|18(l we can write this as 

x c {f c ) = A{f c )[d{f c )-»l(f c )}. (20) 

with [i = g x /(K y g y — K x g x ). Since the calibration line is injected with an amplitude 
that is large compared to the external strain, the residual strain is dominated by the 
control signal at the frequency of the calibration line. So to a good approximation we 
can write the feedback equation for the control signal, Eq. (J5J, as 

Xr(fc) « -X c (J e ). (21) 

Combining Eqs. (121 H and 1)21(1 with the feedback loop equations, © and (jHJ) which 
remain un-changed, we arrive at two independent expressions; one for the product of 
a(t) and (3(t) 

<*(t)fi(t) ™ r)f\i rd[ ! C (f \ (22) 

Gq{Jc) 1 - rd(Jc) 
where = dj and another for a(t), 

A{fc)Co(fc) 1 - r d {f c ) 

where r q = q/fxl. 

Thus, the sensing function gain a(t) as well as its product with the feedback filter 
(3(t) are given by expressions involving the complex functions A, Cq and Go as well as 
complex ratios of the error and digital control signals to the injected excitation. 

The digital feedback filter gain (3(t) is recorded by the data acquisition system 
and so Eqs. I|22|) and 123|) should be regarded as two independent ways of calculating 
the optical gain a(t). 

The ratios of the error and digital control signals to the injected excitation are 
computed as follows. We take a fixed amount of time for each of the three time- 
series q(t), d(t) and l(t), and apply a Hann window. We then perform a complex 
heterodyne at the frequency of the calibration line f c and finally integrate over time. 
This procedure yields complex q(f c ), d(f c ) and l(f c ) which are used to compute the 
ratios r& and r q . 

Since the quantities in Eqs. (|22|l and (|23|l are complex and the measurement 
contains noise from the external strain the optical gain a(t) has a small imaginary 
part. This imaginary part is small only to the degree to which Eq. (|21|l is a good 
approximation. It should be noted however that Eq. lt2"Tl) can be made arbitrarily 
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Figure 3. The magnitude (top left) and phase response (bottom left) of the 
analog sensing function Co(/) (solid curve) and the modified digitised sensing 
function Cq(/) (dotted curve). On the top right the absolute value of the relative 
error between the two magnitude responses and on the bottom right the difference 
in degrees of the two phase responses. 

accurate by either making the amplitude of the injected line arbitrarily large or, since 
the external strain has no component that is correlated with the calibration lines, by 
integrating for an arbitrarily long amount of time. 

The optical gain a(t) is typically calculated on time-scales of a few tens of seconds 
which, given the typical injected amplitudes of the calibration lines, is sufficient to 
reduce measurement noise errors to an acceptable few percent. 

4. Implementation 

4-1. Digitisation of filters and filter performance 

The three filters of the feedback loop can be described in terms of analog zero and 
pole models. This is the starting point of the digitisation procedure. 

The sensing function Co, as we have already discussed, consists essentially of a 
cavity pole at around 90 Hz so the inverse of this filter is a zero at around 90 Hz. 
Unfortunately, digitisation of a zero does not lead to a well behaved digital filter. In 
order to resolve this problem we introduce a pole at high frequencies in the inverted 
sensing function Cq 1 prior to digitisation. To ensure the magnitude and phase 



Making h(t) for LIGO 



8 




10 1 10 2 1 3 1 4 1 1 10 2 1 3 1 4 

Hz Hz 

Figure 4. The magnitude (top left) and phase response (bottom left) of the 
original feedback filter Do (solid curve) and the modified feedback filter D' (dotted 
curve). On the top right the absolute value of the relative error between the two 
magnitude responses and on the bottom right the difference in degrees of the two 
phase responses. 

response of the filter do not change in the frequency range of interest§ the additional 
pole is placed at a sufficiently large value of the frequency. We then digitise the analog 
filter using a bi-linear transform at a sample rate consistent with the presence of the 
additional pole. 

The sensing function also contains a unity gain 8 th order anti-aliasing elliptic 
filter at 7570Hz. This filter introduces frequency dependent phase shifts in the data 
that need to be accounted for. The anti-aliasing filter contains a series of zeros on 
the imaginary axis, which in the inverted filter turn into poles. To ensure that we 
obtain a well-behaved inverse digital filter the zeros need to be moved away from the 
imaginary axis prior to digitisation. 

For the Hanford 4km interferometer in the S2 run we generated a modified digital 
filter Cq with a bi-linear transform at a sample rate 16 times greater than the sample 
rate the of the time series q(t), i.e. 262144 Hz. The zeros of the elliptic filter were 
moved away from the imaginary axis by approximately 160 Hz and a single real pole 
was added at 10 5 Hz. 

In Figure |21 we show the magnitude and phase responses and errors of the analog 
sensing function Cq, as well as our modified digitised sensing function C . The errors 

§ We chose the range to be from 30 Hz to 3000 Hz. This range covers the frequency bands used in 
gravitational wave data analyses performed by the LIGO Scientific Collaboration so far. 
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Figure 5. The magnitude (top left) and phase response (bottom left) of the 
original actuation function A (solid curve) and the modified actuation function 
A' (dotted curve) . On the top right the absolute value of the relative error between 
the two magnitude responses and on the bottom right the difference in degrees of 
the two phase responses. The large errors in the phase beyond a few hundred Hz 
are un-important because the contribution to the external strain h(t) from the 
control signal x c (t) is negligible at these frequencies. 

between 30Hz and 3kHz are less than about 1% in magnitude and 5° in phase. 

The LIGO digital feedback filters Dq are constructed from bi-linear transforma- 
tions at 16384 Hz of a series of analog filters. The main one of these is a double pole 
at DC which, on its own outside a feedback loop, is inherently unstable. In order to 
stabilise it, the analog filter can be modified by shifting the double pole away from 
zero by a small amount prior to digitisation. It should be noted that this procedure 
changes the sign of the feedback filter at DC and introduces errors in the response of 
the filter at low frequencies. The rest of the filters in Do can be used exactly as they 
are implemented in the instrument. 

Figure 21 shows the magnitude and phase responses, as well as the errors, of 
the original digital feedback filter Do and the modified digital feedback filter D' for 
the Hanford 4km interferometer during S2. As expected the errors are larger at low 
frequencies. However, the errors between 30 Hz and 3 kHz are less than about 1% in 
magnitude and 6° in phase. 

The frequency response of the actuation function A is primarily determined 
by the pendulum suspension of the mirrors and thus consists of a pair of complex 
poles at the pendulum frequency f p w 0.75Hz. Additionally, the actuation function 
contains an anti-imaging 4 th order elliptic filter at 7570Hz, a snubber used to suppress 
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resonances in the electronics and a measured lAO/is time-delay modeled using a 4 
order Pade filter. The actuator also contains a series of digital notch filters that remove 
components of the digital control signal at the violin mode frequencies of the mirror 
suspensions. 

All the analog parts of the actuator have been digitised using a bi-linear transform 
at 16384 Hz. We have used the existing digital notch filters as they are implemented 
in the instrument. 

FigureElshows the result of this procedure. Plotted are the magnitude and phase 
response of the original actuation function A and the digitised actuation function used 
in the calibration procedure A' as well as the errors for the Hanford 4km interferometer 
during S2. The errors in the magnitude and phase are less than about 1% and 5° 
respectively at frequencies up to a few hundred Hz. At higher frequencies the errors 
in the phase become larger but this is un-important: the contribution to the external 
strain h(t) from the control signal x c (t) is negligible above a few hundred Hz. 

4-2. Signal processing pipeline 

LIGO data is available in so-called "science segments". At these times the cavities 
are in resonance and the instrument is sufficiently stable to produce data of scientific 
quality. The length of science segments depends mostly on the level of seismic activity 
at the location of each of the sites and can last anywhere between a few tens of seconds 
to tens of thousands of seconds. 

Science runs typically consist of a few hundred science segments. This lends the 
calculation of the external strain for an entire science run amenable to a cost-effective 
Beowulf cluster, which is an ensemble of loosely coupled processors with a simple 
network architecture. In our implementation each processor computes the external 
strain for a science segment. The few hundred processes required can be scheduled 
to run on the cluster with a batch system such as Condor Indeed, the pipeline 
described below was used to compute the external strain for the entire S2 run for 
all three LIGO interferometers using Condor on the 300-node Medusa Cluster [5] at 
UWM. 

The first step of our procedure is to compute a table of values for the optical 
gain a(t) and the digital feedback filter gain (3(t) for an entire science segment. As we 
have described in Section obtaining reliable values of the optical gain requires an 
integration time of a few tens of seconds, so each science segment will have anywhere 
from a few to a few thousand values of a(t). For the Hanford 4km interferometer 
during S2 we chose to compute the optical gain every 60s. 

Science segments can be up to a few tens of thousands of seconds long. Therefore 
we cannot keep all the data for a science segment in the memory of a single node. As 
a result, in our signal processing pipeline, each science segment must be divided into 
a number of shorter parts. In the current implementation of the pipeline the data is 
calibrated in 16 second sub-segments. 

Figure shows the digital processing pipeline used in the computation of the 
external strain h(t) from the error signal, i.e. our implementation of Eq. I|16(l . The 
top and bottom branches show the computations of the residual and control strains 
respectively. 

To avoid problems involving the dynamic range of double precision floating point 
numbers, we first high-pass filter 16s of the error signal q(t). For the Hanford 
4km interferometer during S2 we implemented the high-pass filter using a 10 th order 
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Figure 6. Signal processing pipeline for time domain calibration. Boxes labeled 
x!6 and X 1/16 denote up-sampling and down-sampling respectively. 



Butterworth filter at 40 Hz. This is a conservative choice and the specific frequency 
and order of the filter that is necessary depends on the low frequency character of the 
data. Generally, we expect these choices to depend on the specific science run and 
interferometer. For numerical stability the high-pass filter is divided into second order 
sections. Furthermore, to ensure the phase of the signal is not adversely affected by the 
filtering procedure, the high-pass Butterworth filter is applied to each segment of the 
data once forward and once in reverse. It should be noted that this procedure is acausal 
but since we are not attempting to simulate a physical process this is unimportant. 

Since the data is being calibrated in segments care must be taken at the 
boundaries between consecutive segments. The finite impulse response time of the 
Butterworth filter produces a discontinuity at the beginning and also at the end of 
each segment (because the high-pass filter is applied forward as well as in reverse). To 
ensure continuity across contiguous segments we must allow for extra time at the end 
of each segment. Figure shows a schematic of the procedure we have implemented 
in the case of three consecutive 16s segments. We use one extra second of data at the 
end of each segment which allows sufficient time for the high-pass filter to settle. We 
then then apply the first half of the extra second to the start of the next segment of 
data. The shaded areas in Figure [7| contain data for which the filters have not settled 
and are replaced with data from the previous or next segment. It should be noted 
that this procedure does not remove the few hundred milliseconds of initial ringing of 
the Butterworth filter at the start of each science segment. 

To compute the residual signal x r (t) we first multiply each sample of the high-pass 
filtered q(t) by a suitably interpolated value of the inverse of the optical gain a~ l {t). 
Interpolations are necessary here because the filter gains are typically computed over 
a time-scale of tens of seconds, and are not available for every sample of the error 
signal q(t). We use the table of values of the optical gain created for the entire science 
segment described above to compute an interpolated a _1 (i) for each sample of q(t) 
using a cubic spline. In principle we could also have used band-limited interpolation 
but in practice this makes little difference. 
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Figure 7. Schematic of three consecutive 16s segments of data. To avoid the 
discontinuities produced by the finite impulse response of the Butterworth filters 
in the pipeline, we filter an extra second of data at the end of each 16 s segment. 
We then use the first half of this second at start of the next segment. Shaded 
areas contain data for which the filters have not settled and are replaced with 
data from the previous or next segment. 

The inverse sensing function filter has been digitised at 262144 Hz, and therefore 
the signal must be up-sampled by a factor of 16. The up-sampling procedure is simple: 
15 zeros are added between samples and the resulting time-series is smoothed with a 
low pass filter at a frequency lower than the Nyquist frequency of the original time- 
series q(t). For the Hanford 4km interferometer during S2 we have chosen a 12 th 
order Butterworth filter at 6 kHz which is applied to each data segment once forward 
and once in reverse in second order sections. To ensure continuity across contiguous 
segments we apply a procedure along the lines described above and in Figured adding 
an extra second of data at the end of each 16s segment. 

We then filter the up-sampled signal through the inverse of the modified sensing 
function T^q/^-i. We keep the history of the digital filter across segments to ensure 
continuity. Finally, to produce the residual signal x r (t) sampled at 16384Hz, we 
down-sample by a factor of 16 by picking out one out of every 16 samples. 

In our pipeline the digital control signal is computed by multiplying each sample 
of the high-pass filtered time series q(t) with an interpolated value of the digital filter 
gain f3(t) and then filtering the result by the modified digital feedback filter Td<. Up 
to differences between the original and modified digital feedback filter this produces 
the digital control signal d(t). Note that we could have skipped these steps and instead 
used the digital control signal that is read out from the instrument. 

The low frequency components of d(t) as read out from the instrument and that 
computed by our pipeline are different. This is mostly due to the high-pass filtering 
done on q(t) prior to filtering it through To' . In Figure |5] we show a plot of the 
digital control signal as recorded by the data acquisition system as well as the one 
computed in our pipeline versus time for a segment of data taken by the Hanford 4 
km interferometer during S2. To allow for easier comparison in Figure [S] both the 
digital control signal as recorded by the data acquisition system and the one produced 
by our pipeline have been high-pass filtered by applying 4 second order Butterworth 
filters at 40 Hz consecutively. This removes the low-frequency components and makes 
the agreement between the two signals manifest. 
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Figure 8. Plot of the digital control signal d(t) as recorded by the data 
acquisition system (circles joined by a line) and the digital control signal produced 
by our pipeline (crosses) versus time in seconds for a segment of data taken by 
the Hanford 4 km interferometer during S2. Due to the high-pass filtering done 
on the error signal q(t) prior to filtering it through the modified feedback filter 
T D i the low frequency components of these two signals are different. To allow for 
easier comparison both signals have been high-pass filtered using 4 second order 
Buttcrworth filters at 40 Hz. 



The digital control signal is then filtered through the modified actuation function 
Ta' which produces the control strain x c (t). The histories for both the digital feedback 
filter Te)' o and the actuator Ta> are stored in memory to ensure continuity across 16s 
segments. 

The calibrated strain h(t) is then computed by summing the residual and control 
strains x r {t) and x c {t). 

5. Conclusions 

We have described a method to compute the strain h(t) from the output of the 
LIGO interferometers. This method has been implemented off-line for all three 
interferometers for the S2 run and results have been presented here for the Hanford 4 
km interferometer. The codes that implement the pipeline as well as the digital filters 
are available under the LSC Algorithm Library (LAL ^U]) and the procedure is fully 
automated under Condor 

The pipeline introduces errors associated with the measurement noise involved in 
the calculation of the optical gain a(t) as well as the digitisation and modifications of 
the filters. We expect these errors to be less than about 10% in amplitude and 10° in 
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phase across the detection band. We have tested the pipeline by comparing Fourier 
transforms of time-domain calibrated data with data calibrated in the frequency 
domain and found the differences to be well within our errors. We have also faithfully 
reproduced the digital control signal in the time domain from the error signal using 
our modified servo filter T G > o (see Figure |SJ). 

At this time a third science run (S3) has taken place and a similar off-line 
procedure will be applied to generate strain data for that run. The current plan 
is to place an on-line system at the LIGO sites to generate h(t) in the near future. 
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